Rethinking the Properties of the Quark-Gluon Plasma at T ~ Tc 
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We argue that although at asymptotically high temperatures the QGP in bulk behaves as a gas of 
weakly interacting quasiparticles (modulo long-range magnetism), at temperatures up to few times 
the critical temperature Tc it displays different properties. If the running of the QCD coupling 
constant continues in the Coulomb phase till the screening length scale, it reaches the strong coupling 
treshold as (mo) ~ 1. As a result, the Coulomb phase supports weakly bound Coulombic s-wave 
cc, light quark and even gg states. The existence of shallow bound states dramatically increases 
the quasiparticle rescattering at low energies, reducing viscosity and thereby explaining why heavy 
ion collisions at RHIC exhibit robust collective phenomena. In conformal gauge theories at finite 
temperature the Coulomb binding persists further in the strong coupling regime, as found for A/" = 4 
SUSY YM in the Maldacena regime. 



Soon after the discovery of QCD, it has been found that 
at high temperature T the color charge is screened [1] 
(rather than anti-screened in the vacuum), and the cor- 
responding phase of matter was named the Quark Gluon 
Plasma (QGP). Also, it has been shown, both analyti- 
cally and on the lattice, that at high temperature bulk 
quarks and gluons exist in it as relatively free propa- 
gating quasiparticles, modulo color-magnetic effects that 
are known to be non-perturbative [2]. Although explicit 
perturbative series for thermodynamical quantities were 
found to be badly divergent, it was still hoped that some 
kind of re-summation will make the weak-coupling quasi- 
particle picture work, as the screening would keep the ef- 
fective coupling weak anywhere at T > Tc- The earliest 
suggested QGP signal was a disappearance of familiar 
hadronic peaks - p,uj,(t> mesons - in the dilepton spec- 
tra [3]. Moreover, even small-size deeply-bound cc states, 
ric J/ip, were expected to melt at T > Tc [4,5]. 

In this paper we will show that at T ^ Tc it is not the 
case. The electric effects are still strong, causing pertur- 
bation theory to fail. But before we go into our reasoning, 
let us point out that a motivation comes from two re- 
cent developments, the experimental discovery of robust 
collective effects at RHIC, known as radial and elliptic 
flows, as well as recent advances in lattice simulations at 
T = (1 — 3)Tc. They have led to two related questions we 
are going to address: i. Are there bound states of quasi- 
particles above the QCD phase transition? ii. Why does 
the Quark-Gluon Plasma behave like a good liquid rather 
than a dilute gas of quasiparticles ? 

The former question should not be confused either with 
the non-perturbative magnetic interaction [2] or with the 
issue of "screening masses" for space-like correlators [6]. 
For earlier discussion of real bound states above Tc in the 
framework of high temperature QCD sum rules see [7]. 
Based on chiral dynamics driven by instanton it was 
found that the pi-sigma chiral multiplet should still exist 
as a resonance in the QGP [8] . 

Recent numerical lattice studies have found that, con- 
trary to earlier expectations, charmonium states remain 
bound at least up to T = 1.5 Tc [9]. Moreover, the time- 



direction correlators at finite T found significant devia- 
tions from free behavior at T 3Tc [10], in quantita- 
tive agreement with [8]. Their analysis by the minimal 
entropy method have suggested that pseudoscalars and 
even vector resonances exist above Tc. Below, we will 
show why these states exist, and not just for quarks but 
for gluonic plasmons as well. 

The main idea [11] is that after deconfinement and 
chiral symmetry restoration at Tc, nothing prevents the 
QCD coupling from running to larger values at lower mo- 
mentum scale until it is stopped at the screening mass 
scale. Although at high T the coupling is always small [1], 
at T « Tc it can reach large value ~ 1 . As we show be- 
low, this makes Coulomb binding of the pertinent s-wave 
levels possible, in spite of screening. 

Let us start with charmonium problem first, using the 
non-relativistic Schrodinger equation (and ignoring col- 
lisional broadening). For the standard radial wave func- 
tion x{t) — it has the usual form with the reduced 
mass m = mc/2.We will use throughout this work the 
same in-QGP screened Coulomb potential of the Debye 
form [12] V = -(4as(r)/3r)exp(-M£,r) and use the (T- 
dependent) screening mass as a definition of the length 
unit to have its shape fixed for fixed coupling. In such 
units the equation to be discussed reads 
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with — mcE/M^. The appearance of a bound level 
corresponds to zero binding E = 0, k = 0. If the coupling 
constant does not run and is a constant, all parameters 
of the problem appear in a single combination. Solving 
the equation, one can find the condition for the bound 
state to exist 
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For example, using A/Sas = 0.471 from the vacimm char- 
monium potential and rUc — 1.32 GeV, as Karsch et al 
[5] did long ago, one finds a restriction on the screen- 
ing mass Md < M'j^^'' = 0.37 GeV. Lattice studies of the 
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screening masses have been carried in [13] and found that 
for the near-critical QGP Md = (2.25 ± .25)r. So, for 
T 7^Tc = 170 MeV one expects Md = 0.34 - 0.42 GeV. 
So, the condition (2) is satisfied marginally if at all, and 
these authors concluded that r]c,J/tp may hardly exist 
inside the QGP phase. 

The loophole in this traditional argument is the as- 
sumption that the gauge coupling constant remains 
frozen in the QGP, at the sam,e value it had in the in- 
vacuum charmonium potential. Our main idea is that it 
does not have to be so: in the QGP, after deconfinement 
and chiral symmetry restoration, non-perturbative effects 
are smaller and the charge continues to run to larger val- 
ues. In Fig. 1(a) we compare two such potentials: al- 
though the version with the running coupling leads to 
a smaller potential at smaller distances, it is larger at 
larger r. One cannot of course rely on perturbative for- 
mulae in a strong coupling regime (the alternative will 
be discussed elsewhere [14]), in this work we have de- 
cided to simply freeze when it reaches 1. The cusp in 
the potential reflects on that, and it occurs at about the 
screening length rMjj w 1. 

The modified potential with the nmning coupling leads 
to a more liberal condition for charmonium binding, 
namely Md < 0.62 GeV. This translates into charmo- 
nium remaining bound at 

r<TEe«i.6r,. (3) 

This value agrees well with recent lattice measurements 
[9] mentioned in the introduction, so the modified color 
potential had passed its first test. 

Are there bound states of qq, and gg? For sim- 
plicity, in this work we ignore quark masses, as well as in- 
stantons and the 17(1) axial anomaly, so one can view the 
chirality of light fermions to be permanently conserved 
in the QGP. Chiral symmetry excludes the usual mass 
from being developed, and L, i?-handed quarks propa- 
gate independently. Nevertheless, the propagating quark 
modes in the QGP have some dispersion curves with 
the so called chiral mass, defined as the energy of the 
mode at zero momentum Mg = u){'p = 0), perturbatively 
Mq = gT/\/6 to the lowest order [15]. In weak cou- 
pling there are two fermionic modes: i. with the same 
chirality and helicity the dispersion curve at small p is 
uj = Mg + p/3 + /iMg + ii. with the opposite chi- 
rality and helicity the mode is often called a "plasmino" , 
its dispersion curve has a shallow minimum at p = Q.llgT 
with the energy E^-^i^ = 0.38 .gT shghtly below Mg. The 
latter is related to a new in-matter gluonic mode, the 
plasmon [1]. 

For a general analysis of these modes see [16], where 
in particular both modes should have the same effective 
mass. Lattice data on the dispersion curves for quasi- 
particles remains rather crude, and restricted to essen- 
tially Coulomb gauge [17]. Here we will assume that 
a;2 = p2 ^ ^i^j^ ^j^g following values (at T = 1.5Tc) 

^=3.9 ±0.2 ^=3.4 ±0.3 (4) 



Note that at such T masses arc large and their ratio is 
very different from the weak coupling prediction l/\/6, 
while at T ~ 3T'c and higher they are reduced toward 
perturbative values (which we would not include in this 
work, for simplicity). 

If the dispersion law is known, the effective equation of 
motion suitable for discussion of the bound state problem 
can be obtained by standard substitution of the covariant 
derivatives in the place of momentum and frequency. It 
has the form of the Schrodinger equation if the dispersion 
curve can be parameterized as w = Mg + /2M' + .., 
in general with two different constants M, M' . Both for 
weak coupling and lattice data, such approximation seem 
to be accurate withing several percents. 

Now, to address the issue of binding, we first note that 
if all effective masses grow linearly with T, including the 
screening mass, the explicit T dependence drops out of 
(2) to the exception of the logarithmic dependence in T 
left out in ag ~ \ /\n{T / Kqc d) ■ This is why the re- 
gion of "strongly coupled QGP" turns out to be rela- 
tively substantial. For a qualitative estimate, let us set 
the coupling to its maximum, = 1. The combination 
of constants is (4 x 3.9 T) /(3 x 2.25 T) = 2.3, larger than 
the critical value (2), so one should expect the occurrence 
of (strong) Coulomb bound states. Although the (plas- 
mon) gluon modes are somewhat lighter than quarks in 
(4), their Coulomb interaction has a larger coefficient due 
to a different Casimir operator for the adjoint representa- 
tion, 3 instead of 4/3. As a result the effective combina- 
tion in the potential is Bnigag/MD, which is about twice 
larger than for quarks, and thus the gluons are bound 
even stronger (modulo collisional broadening). 

Solving the equations one can make a quantitative 
analysis, using the same potential as above. We found 
that the highest temperature T at which light quark 
states are Coulomb bound is somehow lower than that 
of charmonium, 

Tgg w 1.45 Tc w 250 MeV , (5) 

while the s-wave gg gluonium states remain bound till 

higher temperatures Tgg « 4Tc. As an example, in 
Fig. 1(b) we show the binding energy of two gluons. 
(In absolute units it reaches about 100 MeV, while it 
is much smaller MeV for quarks.) The wave functions of 
all loosely bound states are similar to that of a deuteron, 
with the usual x ^ e^'^^ behavior outside the potential. 

How many bound states are there? Stating with 
quarks, we note that if one ignores current quark masses 
and instanton effects, the chiralities (left and right) 
are conserved. Furthermore, there are two different 
modes of quarks depending on helicity ("particle" and 
"plasmino"), times antiquarks, times flavors. All in all 
ANj- states, connected by continuity to the pseudoscalar, 
scalar, vector and axial vector nonets in the vacuum (see 
Fig.2). 
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FIG. 1. (a)The combination {Mc/MD)rV{r), versus the 
distance rMo ■ The dashed line corresponds to 4as/3 = 0.471, 
while the solid lines correspond to the running coupling con- 
stant, for T = 1; 1.5; 3; lOTc, from bottom up. The cusps 
occur when as reaches 1. (b) The (modulus of the) binding 
energy of two gluons, versus T /T^. 

Now wc come to gluons: both may have any of 3 polar- 
izations ( with longitudinal "plasmons in matter). Fur- 
thermore (as discussed in our longer paper [14]) there 
can be binary composites even if the total charge is non- 
zero, which provides ~ N"^ states. On top of that, one 
can also think of hybrid (qg) composites, with ^ NfNc 
states. These exotic states however need stronger cou- 
pling and may exist only very close to Tc, if at all. 

QGP kinetics was studied in weak coupling, with a 
rather large mean free path, TImfp ~ l/gHn{l/g) 1. 
Similar pQCD-inspired ideas have led to the somehow 
pessimistic expectation that the Rclativistic Heavy Ion 
Collider (RHIC) project in Brookhavcn National Lab- 
oratory, would produce a firework of multiple jets and 
mini-jets with perturbative corrections for rescattering. 
To the contrary, the very first RHIC run in the sum- 
mer of 2000 has shown spectacular collective motion pat- 
terns, known as radial and elliptic flows. The spectra of 
about 99% of all kinds of secondaries, except the high-pj 
tails, were shown to be described by hydrodynamics [18], 
based on the opposite - zero mean free path - regime. 
Further studies of partonic cascades [19] and viscosity 
corrections [20] have found that the data can only be ex- 
plained if partonic cross sections are increased by huge 
factors, from 2-3 mb for pQCD estimates to about 100 
mb. The equation of state used in [18] is not very different 
from lattice predictions and even simple weak coupling 
expectations. 

A sequence of loosely bound states (black dots in Fig. 2) 
should be very important for the QGP kinetics. The 
scattering cross section may be strongly enhanced by 
low- lying resonances, as known examples from nuclear 
physics (such as low-energy neutron-capture resonances) 
clearly demonstrate. The Wigner cross section of slow 
particles 



has a peak when there is both a shallow level and a vir- 
tual state, as we know from NN scattering in the triplet 
and singlet channels [21]. With the denominator S-|-|Sbl 
ranging from 100 to 10 MeV, the corresponding gg cross 
section is in the range 50-500 mb. These huge cross sec- 
tions make quasi-free propagation of quasiparticles im- 
possible. 

We stress that our analysis is not based on the in- 
vacuum effective potential [22], but on the in- matter 
screened one. The bound states emerge in a screened 
but still strongly coupled Coulomb phase. Still, not 
all of the many-body effects are included in the screen- 
ing. Indeed, the range of distances in the weakly bound 
states is large, and by including the l/N^ factor needed 
for absolute matching of the gluon colors, we may even 
have more than one matching gluon inside the available 
bound-state volume. More detailed many-body studies 
should be made, such as the ones carried for "excitonic 
matter" in semiconductors and insulators. We recall that 
the Coulomb parameter in QED plasmas is e'^/ev where 
e-^ = 1/137 is the usual flne structure constant, e is the 
dielectric constant of the background substance and u is a 
typical particle/hole velocity. It can in principle be tuned 
to be around 1 , resulting in bound excitons (particle- hole 
pairs). The latter occur despite Debye-screening, a sit- 
uation much alike ours. Depending on a number of pa- 
rameters, including the density of excitons and temper- 
ature, the system exhibits various phases, ranging from 
an ideal gas of excitons to a liquid or plasma, or even a 
Bose-condensed gas. On its way from a gas to a liquid, 
clustering with 3- and 4-body states play an important 
role. Although one cannot directly relate these two prob- 
lems (quarks and gluons have Nc and (N^ — 1) colors re- 
spectively, while particles and holes have simply charges 
±e), one may think that in the QGP at T ^ Tc some of 
these phenomena may well be there. 

More generally, strongly coupled many-body problems, 
with large or divergent binary scattering length, are stud- 
ied in at least two other settings: (i) a gas of neutrons, 
with a = 18 fm due to a virtual level; (ii) trapped atomic 
Li^, Li^ in which the scattering length can be tuned. Its 
hydrodynamic properties are of great interest. 

Discussion and Outlook. The main assumption 
made in this work is that at T > To the gauge coupling 
is allowed to run till the screening scale. As a result, we 
found that the QCD phase above Tc supports Coulomb 
bound states. In this work we were mainly interested in 
the location of the point at which the s-wave state disap- 
pear. However, more generally, one may ask how these 
states arc related to know hadrons at T = 0. How it sup- 
posed to happen is shown in Fig. 2 (for 2 massless flavors, 
when the phase transition is second order and all lines 
must be continuous). Very close to Tc the quark (chiral) 
mass is expected to become large, so here mesons are 
rather deeply bound states in this region and both rela- 
tivistic and non-perturbative effects become important: 
we plan to study this region elsewhere. It would also be 
interesting to study whether annihilation of vector states 
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can be seen in the dilepton spectra, and how the survival 
of charmonium above affects its evolution at RHIC. 
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FIG. 2. Schematic dependence of hadronic masses on tem- 
perature T (in units of the critical one Tc), for 2-flavor QCD 
in the chiral limit. The dash-dotted line corresponds to twice 
the (chiral) effective mgiss of a quark. Black dots marked 
s,p,d correspond to the points where the binding vanishes 
for states with orbital momentum Z = 0, 1, 2.... 

In QCD, as soon as the lowest state hits zero there is a 
phase transition. This is however impossible for confor- 
mal gauge theories (CFTs), such as A/'=4 supersymmet- 
ric gauge theory. In CFT the gauge coupling is allowed 
to become supercritical or even large A = g'^Nc » 1. 
Maldacena's AdS/CFT duality [23] has opened a way 
to study this strong coupling limit using classical grav- 
ity. At finite T it was recently actively discussed, for 
Debye screening [25], bulk thermodynamics [24] and ki- 
netics [26] . Although the thermodynamical quantities are 
only modified by an overall factor of 3/4 in comparison 
to the black-body limit, kinetics is changed dramatically. 
In our separate paper [14] we show that in this regime the 
matter is made of very deeply bound binary composites, 
in which the supercritical Coulomb can be balanced by 
centrifugal force, a re-summation of a class of diagrams, 
in vacuum and at finite T. Specific towers of such bound 
states can be considered as a continuation of Fig. 2 to the 
left, toward stronger and stronger coupling. 
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